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ABSTRACT 

Cloud evolution for various metallicities is investigated by three-dimensional 
nested grid simulations, in which the initial ratio of rotational to gravitational 
energy of the host cloud j3 (= lO^-lCT 6 ) and cloud metallicity Z (= 0-Z Q ) are 
parameters. Starting from a central number density of n c = 10 4 cm~ 3 , cloud evo- 
lution for 48 models is calculated until the protostar is formed (n c ~ 10 23 cm -3 ) or 
fragmentation occurs. The fragmentation condition depends both on the initial 
rotational energy and cloud metallicity. Cloud rotation promotes fragmentation, 
while fragmentation tends to be suppressed in clouds with higher metallicity. 
Fragmentation occurs when /3 > 1CT 3 in clouds with solar metallicity [Z = Z Q ), 
while fragmentation occurs when (3 > 10~ 5 in the primordial gas cloud (Z = 0). 
Clouds with lower metallicity have larger probability of fragmentation, which 
indicates that the binary frequency is a decreasing function of cloud metallicity. 
Thus, the binary frequency at the early universe (or lower metallicity environ- 
ment) is higher than at present day (or higher metallicity environment). In 
addition, binary stars born from low-metallicity clouds have shorter orbital peri- 
ods than those from high-metallicity clouds. These trends are explained in terms 
of the thermal history of the collapsing cloud. 

Subject headings: binaries: general — cosmology: theory — early universe — stars: 
formation 



Introduction 



Observations have shown that about 60-80% of field st ars are members of binary or 
multiple systems (e.g., lAbtlll983l ; iDuquennoy &: Mayor Ill99ll ). Hence, a majority of main- 



sequence stars belong to binary or multiple systems. These stars h aye a metallicity equiv- 



alent to that of the sun (i.e., solar metallicity). On the other hand, iLucatello et al.l (120051 ) 
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investigated radial velocities of carbon-enhanced, s-process-rich, very metal-poor (CEMP-s) 
stars that constitute a substantial proportion of the extremely metal-poor (EMP) stars of 
the Galactic halo. In their observations, they found that about 100% of these metal-poor 
stars are in binary (or multiple) systems, and the orbital periods of these stars are shorter 
than those of field stars in the solar neighborhood. Re cently, EMP stars with [Fe/H] < — 5 



(HE0 1 07-5240 and HE132 7-2326) have been observed flChristlieb et al. 11200 it iFrebel et al. 



20051 ). ISuda et al. I (120041 ) showed that the observed pattern of metal abundance for these 
stars can be explained by nucleos ynthesis and mass transfer in a first-generation low-mass 
binary star. iKomiya et al.l (120071 ) investigated the origin of EMP stars in the context of 
stellar evolution and found that the observed EMP stars were exclusively born as secondary 
members of binaries. Thus, these studies imply a high binary frequency in the early universe 
(or low metallicity environment). 

For the present day star formation process, detailed numerical simulations studies of 
fragmentation and binary formation process ha ve been performed by many authors (see 
Bodenheimer et al. Il2000l ; iGoodwin et al. 1 120071 ). The cloud evolution has been calculated 



from n c ~ 10 4 cm 3 to ~ 10 21 cm 3 , where n c is the number density at the cloud center. The 
studies have shown that fragmentatio n occurs only for 10 11 cm -3 <^ n c <^ 10 15 cm -3 (see also 
Bate!ll998l : Iwhitehouse fc Batdl2~00~6L 



On the other hand, in the colla psing primordial cloud (Z = 0), fragmen t ation rarely 



occur s for 10 4 <^ n c <^ 10 16 cm 3 (e.g., lYoshida et al. 



2006 



2002h and frequently occurs for n ^ 10 16 cm 3 ( Saigo et al 



Bromm et al 



2004 



2002 



Abel et al. 



Machida et al.l l2008bh . 



Machida et al.l (j2008bl ) also showed that binary separation in the early universe is narrower 
than that at the present day. Recently, star formation in a c ollaps ing l ow metallicity (Z < 
10~ 4 Zn) gas cloud has been studied by ISmith fc Sigurdssonl (120071 ) and I Clark et al.l (120081 ). 
Clark et al.l (120081 ) found that binary or multiple stellar systems can form even in a low 
metallicity cloud. However, in their studies, since they adopted a sink cell, they calculated 
the cloud evolution only up to n c ~ 1 16 cm" 3 . Fragmentatio n may occur in a later collapsing 



phase (n c ^ 10 cm ), as shown in iMachida et al.l (l2008bf ) 



In this paper, we adopt a barotropic equation of state and calculate the evolutions of 
collapsing clouds with various metallicities [Z = — Z & ), from the formation of dense cores 
(n ~ 10 4 cm -3 ) up to stellar core formation (n ~ 10 23 cm~ 3 ). The calculations indicate that 
the binary frequency increases as cloud metallicity lowers, and binary separation in lower 
metallicity clouds is narrower than in higher metallicity clouds. 
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Model and Numerical Method 



To study the evolution of star-forming cores in a large dynamic range of density and 
spatial scale, a three-dimensional nested grid method is use d and the equations of hydrody- 
namics including self-gravity are solved (see Eq. [l]-[3] of iMachida et al.ll2008bl ). For gas 
pressures in collap sing clouds with differ ent metallicities, barotropic relations that approxi- 
mate the results of lOmukai et al. I (120051) are adopted . Figu re [[] shows these relations (thick 
solid lines) as well as the data from lOmukai et al. I (120051 . thick dotted line) for pressure 
evolution with different metallicities plotted as a function of numb er density . To stres s 
the variations of pressure with density, P/n is plotted (for details see lOmukai et al. II2005I ). 
which is proportional to the gas temperature if the mean molecular weight is constant. 

As the initial state we take a spherical cloud of density 1.4 times higher than the hydro- 
static equilibrium with external pressure (i.e., the so-called critical Bonnor-Ebert sphere). 
The initial central density is taken as n c = 1.4 x 10 4 cm -3 . Each cl oud rotates rigidly (Vt { 



around the z-axis. The initial temperatures, which are derived from lOmukai et al. I (120051 ) 



are different in clouds with different metallicities (see dotted lines of Fig. 1). For example, 
a cloud with Z = (primordial composition) has an initial temperature of 230 K, while a 
cloud with Z = Z Q (solar composition) has 7K as the initial temperature. Since critical 
Bonnor-Ebert spheres are assumed as the initial state, the radii of the initial spheres are dif- 
ferent depending on the initial temperature (or assumed metallicity): the radius of a Z = 
cloud is 5.5 x 10 5 AU, while that of a Z = Z Q cloud is 1.2 x 10 5 AU. We have confirmed 
that these initial differences of radius do not greatly affect the subsequent cloud evolution 
in calculations in which the size and density of the initial cloud are changed. 

The models are characterized by two parameters: the initial rotational energy (3q and the 
cloud metallicity Z. The values used for these parameters are /3q = 10" 1 , 10~ 2 , 10 -3 , 10~ 4 , 
10~ 5 and 10" 6 , and Z = 0, 10" 6 , 10" 5 , 10" 4 , 10" 3 , 10~ 2 , 10" 1 , and 1Z G . Combining these 
two ranges of values, cloud evolution for 48 models is investigated. To induce fragmentation, 
1% of the non-axisymmetric density perturbation of the m = 2 mode (i.e., bar mode) is 
added to the initial cloud. 

To calculate a large spatial scale, the nested gri d method is adopted (for details see 
Matsumoto fc Tomisakal booi IMachida et ~aj71 l2005al fbh . Each level of a rectangular grid 
has the same number of cells (128 x 128 x 8), although the cell width h(l) depends on the 
grid level I. The cell width is reduced by a factor of two for every upper level. The calculation 
is first performed with three grid levels (I = 1, 2, 3). The box size of the coarsest grid / = 1 
is chosen to be 2R C , where R c is the radius of the critical Bonnor-Ebert sphere. A new 
finer grid is generated whenever the minimum local Jeans length Aj becomes smaller than 
8 h(l max ), where h is the cell width. The maximum level of grids is restricted to Z max ^ 30. 
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3. Results 

In each model, starting from a nearly hydrostatic core with central density n c = 
10 4 cm -3 , the evolution of the collapsing cloud is calculated. When fragmentation does 
not occur, the cloud evolution is calculated until the protostar is formed at n c ~ 10 22 cm -3 . 
When fragmentation does occur, the calculation is often stopped after fragmentation (i.e., 
before the protostar formation), because fragments escape from the finest grid. 

Figure [2] shows the final state for each model against the metallicity Z (x-axis) and initial 
rotation energy (3 (y-axis). In the figure, the cloud evolutions are classified into four types: 
fragmentation (red panel border), merger (violet), non- fragmentation (blue), and stable core 
(gray) models. Fragmentation occurs and two or more fragments appear in fragmentation 
and merger models. After fragmentation, fragments survive without merger until the end 
of the calculation in fragmentation models, while fragments merge to form a single core in 
merger models. In non-fragmentation and stable core models, fragmentation does not occur. 
A single protostar is formed in non-fragmentation models, while a long-lived core is formed 
before protostar formation in stable core models. 

The rightmost column in Figure [2] shows the final states for models with solar metal- 
licity (Z = Z Q ). Including radiative effects, the evolution of clouds with solar metallicity 
has been investigated by many authors. Assuming spherical symmetry, the evoluti on from 



the molecular cloud to the stellar core has been calculate d by many authors ( e .g.. lLarson 



1969 



Masunaga fe Inutsuka I l2000h . In three dimensions, IWhitehouse fc Bate! (120061 ) and 



Stamatellos et al.l (120071 ) have calculated stellar core formation from the molecular cloud 
core. These studies have shown that the molecular gas obeys the isothermal equation of 
state with a temperature of ~ 10 K until n c ~ 10 11 cm -3 , then the cloud collapses almost 
adiabatically (10 11 cm -3 <; n c ^ 10 16 cm -3 ; adiabatic phase) and a quasi-static core (here- 
after, first adiabatic core) forms during the adiabatic phase. Thermal evolution for a cloud 
with solar metallicity is also confirmed in Figure [TJ in which the gas temperature increases 
adiabatically after the number density reaches n ~ 10 11 cm" 3 . In models with solar metallic- 



ity, the first adiabatic core surrounded by a shock front is formed at n c = 10 — 10 cm . 

The rightmost column of Figure [2] (models with Z = Z & ) shows that fragmentation 
occurs when the initial cloud has (3q > 10" 3 . In these clouds, fragmentation occurs only 
after the first adiabatic core formation: fragmentation occurred at n c = 6.8 x 10 11 cm -3 
(fa = lO" 1 ), n c = 2.4 x 10 12 cm~ 3 (fa = 10" 2 ), and n c = 1.7 x 10 13 cm" 3 (fa = 10~ 3 ), 
respectively. In these fragmentation models, fragments survived without merger until the end 
of the calculation for models with fa = 10 _1 and 10~ 2 , while a single stable adiabatic core was 
formed after merger in the model with f3 = 10~ 3 . The frag mentation condition {(3n > 10 ~ 6 ' 



for a cloud with solar metallicity is consistent with that of iMatsumoto fc Hanawal (j2003bl ). 
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Matsumoto fc Hanawal (l2003bl ) calculated the evolution of a cloud with solar metallicity in 



a large parameter space and found that fragmentation occurs when the initial cloud has 
/?o > 2.2 x 10~ 3 . In addition, many st udies have shown that fragmentation occurs only after 
the first adiabatic core formation (see lGoodwin et al. 1120071 ). 



For clouds with Z = Z®, fragmentation occurs in models with (3q > 10 -3 after the 
first adiabatic core formation, while stable first adiabatic cores are formed in models with 
(3 < 1CT 4 . In models with /3 < 10~ 4 , although the cloud evolution was calculated for 
a sufficiently long time, the first adiabatic core did not collapse to reach the protostellar 



density 



10 cm 3 ). When a non-rotating cloud is adopted as the initial state, the 



first adiabatic core incr eases its mass with time by gas accretion, and can collapse again 
(the second collapse, see iMasunaga &: Inutsuka 1 120001 ) to form a protostar for a short dura- 
tion. On the other hand, a long-lived first a diabatic core is frequently formed in a rotating 
cloud, as shown in ISaigo fc Tomisakal (120061 ). This long-lived core can collapse further af- 
ter the non-axisymme tric perturbation grows sufficiently owing to the bar mode instability 
(IDurisen et al. 1 119861 ) , beca use the angular momentum of the core is transferred by the 
non-axisymmetric structure (]Batdll998l : ISaigo fc Tomisakall2006l ). To investigate further the 
evolution for models with (3 < 10~ 4 , in different models of Figure [21 a large initial ampli- 
tude of the non-axisymmetric perturbation is adopted. In these calculations, after the first 
adiabatic core formation, the non-axisymmetric pattern grows and the cloud collapses (the 



second collapse) to reach the stellar density 



Tic 



10 



22 



cm 3 ). However, these models d o not 
£ 10 14 ' 



show fragmentation even in the later evolution phase (n c ^ 10 cm ), as shown in [Bate 
(119981 ) . Thus, it is expected that fragmentation does not occur in stable core models even 
after a calculation for an extended time period. This is because the central region transforms 
from a disk-like to a spherical configuration due to the large thermal energy after the first 
adiabatic core formation; fragmentation can easily occur in a thin disk-like configuration. In 
addition, the angular momentum is effectively transferred from the high-density region by 
the non-axisymmetric pattern. 

The leftmost column of Figure [2J shows final states for models with Z = (primordial 
cloud). In the primordial cloud, fragmentation oc curs in clouds with (3q > 10~ 6 , which is con- 
sistent with the results of iMachida et al.l (l2008bl ). These panels also show th at the fragmen- 
tation scale (i.e., distance between fragments) increases with /3 . As shown in lMachida et al. 



d2008al lbh. the fragmentation scale is comparable to the Jeans length at the fragmentation 
epoch. Since fragmentation occurs in the earlier epoch (or lower density) for models with 
larger (3 and the Jeans length shortens as the density increases, fragments have a larger 
separation in clouds with larger j3 . 



The first adiabatic core is formed at n c 



10 11 -10 14 cm 3 in clouds with solar metallicity, 
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while the protostar is formed directly without the first adiabatic core formation in the pri- 
mordial cloud. As shown in Figure [U since the temperatu re in the primordial clou d (Z = 0) 



gradually increases with 7 ~ 1.1 of the polytropic index (jOmukai fc Nishi 1 119981 ) in a wide 
density range of 10 3 cm~ 3 < n c < 10 18 cm -3 , the (first) adiabatic core is not formed before 
the protostar formation (n c <^ 10 21 cm -3 ). Thus, fragmentation can occur for n c <; 10 21 cm -3 
in the primordial cloud, while fragmentation occurs only for n c <^ 10 14 cm" 3 in clouds with 
solar metallicity. In addition, in the primordial cloud, since the central region can be spun 
unti l the protostar is form ed (n c ^ 10 21 cm -3 ) due to the absence of the first adiabatic core 



sec 



Machida et al.l 120071 ). fragmentation occurs even in clouds with smaller /Jo- 
in Figure [2l the border between fragmentation and non-fragmentation models is shown 
by a gray-white dotted line, which indicates that fragmentation tends to occur in clouds with 
lower metallicity. For example, fragmentation occurs in models with (3q > 10 -4 for clouds 
with Z = 1O~ 3 Z and 10~ 4 Z Q , while fragmentation occurs for (3q > 10~ 3 for clouds with 
Z = IO^Zq and 1O~ 2 Z . 



The fragmentation epoch is closely related to the first adiabatic core formation epoch. 
The first adiabatic core is formed in clouds with Z > 10~ 5 Z Q , while the first adiabatic core is 
not formed in clouds with Z < 10~ 5 Z & . As shown in Figured], the gas temperature increases 
keeping 7 ~ 1.1 for clouds with Z = and 10~ 6 Z & , while the gas temperature increases 
adiabatically (7 ~ 1.4) after the central density reaches n c ~ 10 n -10 16 cm -3 in clouds with 
Z < 1O _5 Z . In clouds with Z < 1O~ 5 Z , the epoch when the cloud collapses adiabatically 
depends on the cloud metallicity. The first adiabatic core is formed at an earlier epoch (or 
lower density) in clouds with lower metallicity. For example, the first adiabatic core is formed 
at n c ~ 5 x 10 11 cm -3 in clouds with (Z, /3q) = (10 _1 , 10~ 6 ), while the first adiabatic core is 
formed at n c ~ 6 x 10 14 cm" 3 in clouds with (Z, (3q) = (10~ 5 , 10 -6 ). When the first adiabatic 
core is formed at a later evolution phase (or higher density), the cloud is spun over a long 
duration and forms a thin rotating disk even in clouds with small (3q. Then, fragmentation 
occurs in the thin disk. 



4. Discussion 



Although we cannot observe the primordial gas cloud i n the early u niverse, we have 
observed many molecular clouds in the solar neighborhood. ICasellil (120021 ) observed about 
60 molecular cloud cores and found them to ha ye rotational energies in the range j3o = 10~ 4 - 



0.07, with a typical value of (3q = 0.02 (see also lGoodman et al. Ill993l ). Figure [2] shows that 



fragmentation occurs in clouds with solar metallicity (Z = Z Q ) when the rotational energy 
exceeds (3q > 10~ 3 , which is smaller than the typical value of the observations. Thus, it is 
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expected that fragmentation frequently occurs in clouds with solar metallicity. 



Observations have shown a high binary frequency at the present day (see lMathieulll994j ) . 
However, we cannot determine the binary frequency in the early universe (or binary frequency 
in a lower metal environment) from observations. Cosmological simulations show that the 
first collapsed objects formed at z ~ 10 with n c ~ 10 3 — 10 4 cm -3 have rot ational energies 



comparable to or slightly l arger than the typical value of molecular clouds (IBromm et al. 



20021 ; lYoshida et al. 1 120061 ) . Thus, if the distribution of cloud rotational energy does not 
strongly depend on metal abundance, the binary frequency will increase with decreasing 
cloud metallicity. Even if a typical low-metallicity cloud has a smaller rotational energy of 
flo <C 0.02, the binary frequency may be high, because fragmentation can occur even in these 
clouds; fragmentation occurs in clouds with Z < 1O _4 Z when (3 > 10~ 4 , which is much 
smaller than the typical value of molecular clouds (/3 = 0.02). Thus, it is expected that the 
binary frequency in low-metallicity clouds is higher than in solar-neighborhood clouds. 

The fragmentation scale or separation between fragments decreases with decreasing 
cloud metallicity. For example, at the fragmentation epoch, the separation between frag- 
ments for model (Z, (3 ) = (Z Q , 0.1) is ~ 45 AU, while the separation for (Z, /3 ) = (0, 0.1) 
is 0.2 AU (see Fig. [2]). Since the first adiabatic core is formed at a higher density in clouds 
with lower metallicity, the fragmentation scale decreases with decreasing metallicity, i.e. the 
fragmentation scale (or Jeans scale) shortens as the cloud density increases. Thus, at the 
moment of birth, the binary orbital period in low-metallicity clouds is shorter than that of 
solar-neighborhood clouds. 

In the models shown in Figure [21 1% of the non-axisymmetric perturbation is added to 
the initial state. When different initial amplitudes of the non-axisymmetric perturbation are 
added, the cloud evolution may change. For example, in the fragmentation region of Figure^ 
although fragments merge to form a single core only in model (Z, (3q) = (10~ 3 , 0.1), fragments 
do not merge in this model when different amplitud es of the non-axisymmetric perturbation 



(5, 10, 20, 30%) are added. iGoodwin et al. I (120071 ) pointed out that statistical studies are 



needed to understand fragmentation, because the cloud evolution in the high-density region 
is highly chaotic. In subsequent papers, we will statistically investigate the cloud evolution 
for different initial cloud shapes in a large parameter space. 



I thank K. Omukai for supplying the data of thermal evolution for the primordial 
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out on VPP5000 at the Center for Computational Astrophysics of the National Astronomical 
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n (cm' 3 ) 



Fig. 1. — Thermal evolution of collapsing clouds against gas number density with different 
metallicities (Z = 0, 10" 6 , 10" 5 , 10" 4 , 10" 3 , 10" 2 , 10 -1 , 1Z & ). To stress the variation of 
pressure with density, P/n is plotted, which is proportional to the gas te mperature if the 



mean molecular weight is constant. The thick solid lines show data from lOmukai et al. 
( 120051 ). The thick dotted lines are fits by our numerical simulation. 
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Fig. 2.— Outcomes of the collapse of prestellar clouds against different metallicities Z 
and initial rotation energies /3q. The density distribution (color-scale) around the center on 
the equatorial plane is shown in each panel. The grid level is displayed at the upper left 
corner of each panel and the grid scale is denoted at each lower right corner. The panel 
border color indicates the following: red: fragmentation model, violet: merger model, blue: 
non-fragmentation model, and gray: stable core model. The gray-white dotted line shows 
the border between fragmentation and non-fragmentation models. 



